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In this work we study the properties of the atomic entanglement in the eigenstates spectrum of the 
inhomogeneous Tavis-Cummings Model. The inhomogeneity is present in the coupling among the 
atoms with quantum electromagnetic field. We calculate analytical expressions for the concurrence 
and we found that this exhibits a strong dependence on the inhomogeneity. 
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I. INTRODUCTION 

Quantum correlations have played a central role from 
the very beginning of Quantum Mechanics, and currently 
the concept of Entanglement has become a key resource 
in the research on quantum information and quantum 
computation The availability of entangled quantum 
states and their characterization are among the most 
important questions in quantum information. In this 
sense, there have been a considerable number of theo- 
retical works that have allowed for characterizing quan- 
tum entanglement. From a mathematical point of view, 
a quantum state of a pair of particles is separable if this 
can be written as a convex sum of product states Q, 

P = J2i Pi\ a i) ( a i\ ® \h) (h\, where |a 4 ) and corre- 
spond to different particles in the pair and pi = 1. Ac- 
tually, separability criteria 3, 4] provides us with neces- 
sary and sufficient conditions to state if a bipartite quan- 
tum state is separable or not. However, these criteria do 
not allow us to quantify exactly the amount of Entan- 
glement of a given system. Bennet et.al. [f| have char- 
acterized the necessary and sufficient channel resources 
to transmit quantum states. This leads to a measure of 
entanglement that is called Entanglement of Formation. 
In other important work, Wootters @, 0] has shown that 
the Entanglement of formation of an arbitrary state of 
two qubits is defined by an exactly calculable quantity 
called concurrence. 

On the other hand, recent advances in the manipula- 
tion of collections of atoms have led to the possibility of 
considering inhomogeneous coupling, for example, in sys- 
tems of atomic clouds. Some contributions have emerged 
inspired by this fact. For instance, decoherence of col- 
lective atomic states due to the inhomogeneous coupling 
between the atoms and external fields has been studied 
by Sun et.al. Q. Also, Kuzmich et.al. @ has pointed 
out the nonsymmetric character of the entanglement of 
multi-atom quantum states due the inhomogeneity. Re- 
cently, decoherence in process of quantum information 
storage in atomic clouds has been studied [lfj, where 
the phenomenon of electromagnetically induced trans- 
parency (EIT) is present. 

Atomic cloud physics could be conveniently described 
by the Dicke model [ll| , when considering the interaction 
of atoms with light in free space. The Tavis-Cummings 



(TC) model [12| , is suitable when the coupling takes place 
inside a cavity. In most applications of the TC model, a 
constant coupling between the atoms and the radiation 
field is assumed. This simplification is, of course, essen- 
tial when an analytical description of the coupled system 
is required [H, M, OJ, OJ] , at least when a small number 
of atoms is involved. For the case of many atoms, the dy- 
namics can still profit from the associated SU(2) group 
structure and numerical solutions can be found. How- 
ever, the situation is drastically different when we con- 
sider the more realistic inhomogeneous coupling. In this 
case, there is no possibility to access the Hilbert space in 
a simple manner, because all angular momenta represen- 
tations are mixed along the dynamics and no analytical 
approach is known. 

The entanglement properties of the ground state for 
the Dicke model has recently been studied by authors in 
Ref. [TtI EH- In such case, when the coupling between 
atoms and field is homogeneous, the concurrence between 
two atoms is independent of the pair considered due to 
the symmetry of the Dicke states. A similar situation 
should arise in the case of the Tavis Cummings model. 
From the point of view of the model properties, it would 
be interesting to study the quantum correlations between 
atoms under more general conditions, as is for the case 
of inhomogeneous coupling. An interesting question to 
answer would be, how does the spatial profile of the cou- 
pling reflect on entanglement between atoms?. 

In this work we aim to describe the entanglement prop- 
erties in the spectrum eigenstates of the TC model [l2[ 
by considering inhomogeneous coupling of the atoms to 
the quantum electromagnetic field. Specifically we study 
the bipartite atomic concurrence by tracing out N — 2 
particles. This work is organized as follows: in Sec. II 
we present the inhomogeneous Tavis-Cummings (ITC) 
model and the properties of the Hilbert space are shown. 
In Sec. Ill, we study the properties of the eigenstates of 
the model, the bipartite concurrence between atoms is 
analyzed for different number of excitations in the sys- 
tem. In Sec. IV, we present our concluding remarks. 



II. THE MODEL 

The Hamiltonian describing the inhomogeneous inter- 
action of N atoms with a single mode of a quantum elec- 



tromagnetic field in the on resonance regime is given by 
H = H + Hj where 



H = a^a + S z , 
Hj = aJ + + cv J_, 



(1) 



besides we have taken H = 1, S z = X^=i °z/2 an d 
J± = Y^,jLi K j a ~L are the atomic collective operators with 
Kj = cjj/w. The number of excitation operator a^a + S z 
is an integral of motion for the system [T^ |. In addi- 
tion Kj is the inhomogeneous coupling of atoms to the 
field. In this work, we will consider Kj = Ksin(7TXj / L). 
We notice that due to the inhomogeneous coupling, the 
collective atomic operators do not satisfy a closed Lie al- 
gebra su(2), leading to a larger total Hilbert space, so 
that, the description in terms of symmetric Dicke states 
is not allowed anymore. 

An effective approach to study the dynamics of systems 
with inhomogeneous coupling can be developed [19j | . The 
basic idea of this approach is to follow the Hilbert space 
that the coupled system will visit along the evolution and 
to implement a truncation criteria based on a probabilis- 
tic argument. In the homogeneous case the eigenstates of 
H are built as superpositions of eigenstates of Ho , that 
is, states with a fixed number of excitations. Thus, the 
ground state corresponds to the state with all atoms in 
the lower state. The first excited state corresponds to a 
superposition of one photon state with all atoms in the 
ground and the state with zero photons and the sym- 
metric state with one excitation. In the inhomogeneous 
case we have to find which are the states with a fixed 
number of excitations. For example, let us consider the 
state |n)| 0) where |n) denotes the state with n photons 
and |0) denotes the collection of N atoms in the ground 
state \gxg2...gpf). Let us see now how the state |n) |0) 
is coupled to other states under the interaction Hamilto- 
nian ([1]) which conserves the number of excitations. We 
see that the state |0)|0) is not coupled by the interaction 
Hamiltonian to other states of the global system, that is, 



|0)|o) — |o>|6>. 



(2) 



Now, if we have one excitation in the electromagnetic 
field |1)|0) this state is coupled through the term J + a so 
that 



TV 



J+o|l)|0) =10)^^11,), 



(3) 



where represents the state where the i-th atom is 
excited. We define the normalized state 



ID = 



l 



TV 



(4) 



such that the matrix element of the hamiltonian between 
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FIG. 1: Levels scheme that involve the Hilbert space. Here is 
possible to see the Hilbert space's growth due to the Inhomo- 
geneous coupling. 



|1)|0) and |0)|1) states is given by 



(0|(I|H|1)|0) 



A 



N 

E< 

i=l 



(5) 



As we can see, the number of excitations is conserved. 
The resulting state is coupled to the state |1)|0) such 
that, starting from this root the system evolves in a closed 
two dimensional subspace {|1)|0),|0)|1)}. If we increase 
the number of excitations k in the system, for exam- 
ple k = 2, one would expect that the Hilbert space will 
be spanned by the vectors {|2)|0),|1)| 1),|0)]2)}, where 

I 2)=^2E i<j «i«j|2tf) and N 2 = ^4£~^ k\k). How- 
ever, it can be shown that the state |0)|2) is cou- 
pled through J_ to the state |1) |$i) , where |$i) = 



which is different from so 



n 2 

we can say that |$i) has a component along |1) and a 
component along a state |l p ) perpendicular to the state 
|1) [HJ. Thus, if the number of excitations keep grow- 
ing, then it is possible to generate a sequence of collective 
atomic states and its respectively perpendicular states as 
shown by Fig. [U where the coupling among atomic col- 
lective states are shown. 

It is possible to show that the states |n) (first row in 
Fig. [U can be obtained from the repeated application of 
the operator aJ+ over the state \n) |0) 



(a J+) n |n) |0) = VW!|0) £ 



^l<^2<■■<^n 



(6) 

this implies that the collective atomic states |n) are given 
by 



(7) 



ll<l2<--<tn 



with N% = ( n -) 2 K i 1 K i 2 -- K i n - O n tn e other 

il<i2<--<in 

hand, the sequence of states |n p ) in the second row of 
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FIG. 2: Contribution (%) of second row of Fig. Q] to the 
energy as a function of the number of atoms N. Ai5' fc - ) = 
W0(\(E[ k) - E ( 2 k) )/E[ k) \), where E (k) 2 are the energies for k 
excitations using 1 or 2 rows of Fig. ff] 



Fig. Q] are given by 



1 



(<!> n \<t> n )-\{n\$ n )\ 2 
where the states |<&„) are written as 



$„) - (n | $„) |n)] , 
(8) 



(n-1)! £ 



n. 



>>> ( 9 ) 



and 



= E 4, 



E < 



N 



4=1 



(10) 



An accurate description of the model should consider all 
the states we have found associated with the inhomo- 
geneity. However, based on numerical calculations we 
have obtained that the spectrum of the Hamiltonian is 
not sensitive to perpendicular states further than the first 
row. This is shown in Fig. [2l where the effects of the sec- 
ond row on the energies obtained from the first row are 
quantified. In this figure we observe that the contribution 
of the second row in Fig. [1] is less than 1%. Although the 
contribution of the second row in Fig. [1] is higher when 
the number of excitations increases, decreases when the 
number of atoms is increased, as can be seen in Fig. [5] 
In other words, the contribution of the second row to the 
energy, will increase only when k/N increases, where k is 
the number of excitations. Thus, we can conclude from 
Fig. [2] that only the first row of Fig. Q] is needed in or- 
der to characterize the spectrum of the ITC model for 



the parameters used in this work. On the other hand, 
it is important to note that this technique of taking into 
account only a few rows in Fig. Q] is also valid when the 
quantum dynamics of the system is studied. In that case 
the approximation is valid within a given window of time. 
The size of this window of time will depend on both TV 
and k El. 



III. EIGENSTATES SPECTRUM PROPERTIES 

In the previous section we have discussed the necessary 
ingredients to deal with the ITC model. Particulary, we 
have shown that only a part of the Hilbert space associ- 
ated to the system is necessary for a suitable description. 
In this section we study the properties of the eigenstates 
of the system, by presenting some analytical calculations 
that allow us to characterize the behavior of the entan- 
glement between a pair of particles. The spectrum and 
bipartite atomic Concurrence Q , can be analytically cal- 
culated in the inhomogeneous case for 0, 1 and 2 excita- 
tions respectively. The energies for these cases are given 
by: 



N 
2 



Ni, 



(11) 



£(2) = 2 - f - J2N? + (Ay/Vr) 2 , 



where (Ai) 2 = J^Zi K ? and (^a) 2 = 4 J2?<j rftf, cor- 
respond to the normalization of the states with one and 
two excitations in Eq. ([7]) respectively. These energies 
are associated with the eigenstates 



^(0) 



= |0)|0), 

= -L{|i)|o)-|o)|T)}, 

= a|2)|6)--^|l)|l)+/5|0)|2) 



(12) 



where 



V2 



Nr 



(3 = 



y/4Nf+2N% ' 
N 2 



(13) 



From the expressions for eigenstates in Eq. ([12]) we are 
able to calculate an analytical expression for the bipartite 
concurrence of a pair of atoms i, j by tracing out respect 
to the other A — 2 atoms. For this cases the concurrence 
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FIG. 3: Bipartite concurrence Cij between the first and the 
j'-th atom for k = 1 excitation for different number of atoms 
N. 



FIG. 4: Bipartite concurrence Cijv between the first and the 
last atom for k = 1,2, ...,6 excitations as a function of the 
number of atoms N. 



is given by: 



a 



(0) 



a 



(i) 



\ K i K j I 



a 



(2) 



--2\k. iKj \[^ + ^M? 



(14) 



iv 2 



l 1/2 ], 



With Ml = Ylkjtij K h and M 2 = T,L<njti,j ^n- 

It is important to note from these expressions for the 
concurrence that while in the homogeneous TC model 
the bipartite concurrence of a pair of atoms is indepen- 
dent of the chosen pair, the apparition of inhomogeneity 
in the coupling between the atoms to the quantum field 
causes the bipartite concurrence to be dependent on the 
pair of atoms that we choose. As is apparent from these 
equations, the entanglement between a pair of atoms is 
proportional to the coupling constants of atoms to the 
field. Thus a pair of atoms located in a region of strong 
coupling to the field will be more entangled. As we will 
see these dependence of entanglement on the coupling 
constant appears considering more excitations. 

Fig. [3] shows the bipartite concurrence C\j, i.e., the 
concurrence between the first and the j-th atom for dif- 
ferent total numbers of atoms N. All atoms are equally 
spaced for each case. We observe in this picture the ef- 
fect of the inhomogeneous coupling of atoms with the 
field. The concurrence exhibit a reminiscence of the spa- 
tial profile of the single mode electromagnetic field. The 
maximum of entanglement happens for the strongest cou- 
pling of the j-th atom, which occur for this atom located 
at the center the cavity. Comparing the entanglement 
of the first atom with the atom located around the cen- 
ter of the cavity, we realize that as the number of atoms 
N inside the cavity is increased, entanglement decreases. 



i Pi 



This is in agreement with that found in the homogeneous 
case. 

In the most general case, the state with k excitations 
can be written as 



*«)=5> s \s) A \k-s) F , 



(15) 



s=0 



where |s)^ is the collective atomic state with s excitations 
defined in Eq. ((TJ) and A s is a coefficient arising from the 
eigenvalue problem. Tracing out the quantum field, the 
reduced density matrix for the atomic subsystem is given 
by 



(16) 



Now, tracing out N — 2 atoms, the bipartite concurrence 
between atoms i,j will be given by 

k \A> 2 



s=0 s 
s,s'=0 s s ' 

where we have defined 



(17) 



Mr, 



ni<ri2<..<n p 



2 



(18) 



E ("0 «ni«na 

ni<n2<-.<n p =£i,j 



As we can see in the general expression for the concur- 
rence flTTl) for k excitations and N atoms, the entangle- 
ment between two of these atoms still depends directly 



5 



o ° ° o 
o o 




o k=1 
□ k=2 


o o 




x k=3 


o o 




• k=4 






k=5 





o 


* k=6 


o 




o 


n D ° D ° r, 
□ v v □ 
□ x x x 




o 


> p X $ * * £ 


□ 

* 


□ ° 


}* 


ft ° 

5 



FIG. 5: Bipartite concurrence Ci,j between the first and the 
j-th atom for k = 1, 2, 6 excitations with iV = 20 atoms. 

on the coupling constant. As shown by this result this 
is one important feature of entanglement that the TC 
model exhibit for the inhomogeneous case. This makes 
the difference with the the situation for the homogeneous 
case where only depends on the total number of atoms N 
for a given number of excitations k [l7j |. Fig.Q]shows the 
concurrence between the first atom and the last atom as a 
function of the number of atoms inside the cavity. In this 
picture we can observe how the concurrence decreased 



when the number of total excitations k in the system in- 
creases. However, when N ^> 1 the concurrence becomes 
less dependent on the number of excitations tending to a 
fixed value. Fig. [5J as Fig. [3] shows the concurrence C±j 
between the first and the j-th atom in the cavity. The 
effects of the inhomogeneous coupling and the loss of en- 
tanglement when the number of excitation is increased 
are also clear in this figure. 



IV. SUMMARY 

In the present work we have found that although the 
symmetry SU(2) no longer exits in the ITC model and 
consequently a treatment in terms of symmetric Dickc 
states is not feasible, there are still some features that 
can be analyzed. In particular, the bipartite concurrence 
of the eigenstates of the model can be obtained exhibiting 
an explicit dependence on the inhomogeneity in addition 
to the dependence on the number of atoms and the num- 
ber of excitations. 
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